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ABSTRACT 


Gasdynamic  properties  which  include  the  effects  of  dissociation 
and  intermolecular  forces  (van  der  Waals)  are  presented  for  incident 
and  reflected  shock  waves  in  equilibrium,  imperfect  nitrogen.  Charts 
are  presented  for  incident  shock  Mach  numbers  in  the  range  from  1 
to  10  into  (ideal)  nitrogen  at  a  temperature  of  300°K  and  pressures  in 
the  range  from  10"^  to  10^  cm  Iig.  The  temperature  and  density  in 
any  region  do  not  exceed,  respectively,  15,  000°K  and  1000  amagats. 

In  addition  to  the  usual  incident  and  reflected  shock  properties,  stag¬ 
nation  conditions  upstream  and  downstream  and  conditions  imme¬ 
diately  downstream  of  a  standing  shock  wave  are  also  presented.  At 
pressures  above  one  atmosphere  in  the  undisturbed  gas,  the  effects 
of  the  intermolecular  forces  on  the  gasdynamic  properties  are  demon¬ 
strated. 
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NOMENCLATURE 

a  Speed  of  sound,  ft /sec 

h  Enthalpy,  ft^/sec^ 

M  Mach  number 

Ms  Incident  shock  Mach  number.  Us/a, 

P  Pressure,  atm 

q  Sphere  stagnation  heat-transfer  rate,  Btu/ft^-sec 

R  Gas  constant,  ft^/sec^  °K 
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rn  Sphere  (nose)  radius,  in. 

S  Entropy,  ft2/sec2  °K 

T  Temperature,  °K 

Ur  Reflected  shock  wave  velocity,  ft/sec 

Us  Incident  shock  velocity,  ft/sec 

u  Gas  velocity,  ft /sec 

ff  Molecular  weight 

Z  Compressibility  factor 

y  Ratio  of  specific  heats 

p  Density,  amagats 

p  Viscosity  coefficient,  dimensionless 

SUPERSCRIPTS 

(n)  Iteration  index 

*  Dimensional  pressure  in  lbf/ft2  and  dimensional  density  in 

lbf-sec2/ft^ 

SUBSCRIPTS 

1  Conditions  in  the  quiescent  gas 

2  Conditions  behind  the  incident  shock 

5  Conditions  behind  the  reflected  shock 

a  Conditions  at  one  atmosphere  pressure  and  273.  15°K 

w  Conditions  at  the  wall  (300°K) 

DOUBLE  SUBSCRIPTS 

ij  Unless  otherwise  noted,  a  ratio,  e.  g. ,  p-J  p^ 

20  Stagnation  conditions  upstream  of  a  standing  shock  wave 
20'  Stagnation  conditions  downstream  of  a  standing  shock  wave 
2s  Conditions  downstream  of  a  standing  shock  wave 
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USEFUL  CONSTANTS 

0a  ~  1105.55  ft/sec  =  0.337006  mm/gsec  or  km/sec 
pa  =  1  atm  =  14.696  Ibf/in.3  =  2116.224  lbf/ft3  =  1.013236 
=  760  mm  Hg  =  1.013246  x  10*  newtons/m3 
R  =  3196.793  ft '/sec3  °K  =  1.98717  cal/mole  °K  -  296.992 
=  7.0943286  x  10“  cal/gm  °K  or  Btu/lbm  °R 
Ta  =  273.15  °K  =  491.67  °R 
W  =  28.0134  gm/mole 

Pa  =  2.423516  x  10  3  Ibf-eecVft*  or  slugs/ftS  =  7.7974207 
=  1.2490294  kg/m’  =  1.2490294  x  10-5  gm/cm* 


REFERENCE  QUANTITIES 


T,  =  300  °K  =  540 °R 


aj  =  1158.729  ft/sec  =  0.35318015  mm/^sec  or  km/sec 
h4  =  yRTl/(y“l)  =■  3.3566326  x  10*  ft2/sec2 
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SECTION  I 
INTRODUCTION 


Lack  of  adequate  thermodynamic  properties  of  gases  at  high  pres¬ 
sure  and  temperatures  up  to  15,  000°K  has  prompted  the  AEDC  to  support 
work  on  these  properties  at  the  National  Bureau  of  Standards  (NBS)  for 
several  years  {Refs.  1  to  4).  Also,  the  specific  heat  and  speed  of  sound 
data  at  these  conditions  have  been  computed  by  Landis  and  Nilson(Ref.  5) 
and  Lewis  and  Neel  (Refs.  6  to  8).  The  more  recent  NBS  data  were 
interpolated  at  constant  entropy  and  at  constant  pressure  by  Neel  and 
Lewis  in  Refs.  8  to  12  (see  also  Ref.  13). 

The  recent  data  of  NBS  and  VKF  noted  above  were  used  by  Lewis 
and  Burgess  to  compute  normal  shock  wave  properties  in  imperfect*  air 
(Ref.  14)  and  nitrogen  (Ref.  15)  ir.  the  range  Mg  =  6  to  30  over  the  pres¬ 
sure  range  =  10-4  to  1000  cm  llg  (see  also  Refs.  16  and  17).  These  im¬ 
perfect  gas  results  were  compared  with  the  perfect  gas  results  of  Feld¬ 
man  (Ref.  18)  for  air  and  Bernstein  for  nitrogen  (Ref.  19).  The  com¬ 
parisons  shov/ed  non -negligible  effects  on  such  quantites  as  Z,  p ,  and  a 
in  the  stagnation  and  reflected  regions. 

The  need  for  data  below  1500C'K  has  prompted  the  VKF  to  extend  the 
NBS  results  to  lower  temperatures  (see  Refs.  20  and  21).  Recently  the 
data  for  nitrogen  have  been  extended  to  higher  pressures  by  Grabau  and 
Brahinsky  (Ref.  22)  using  an  extrapolation  procedure  and  requiring  in¬ 
ternal  consistency  and  satisfaction  of  the  laws  of  thermodynamics. 

The  recent  NBS  data  (Ref.  3),  the  low  pressure  data  of  Humphrey 
and  Neel  (Ref.  20),  and  some  recent  unpublished  data  of  Grabau  and 
Brahinsky  using  the  methods  described  in  Ref.  22  were  used  by  Lewis 
and  Burgess  (Ref.  23)  to  compute  the  normal  shock  wave  properties  in 
imperfect  air  in  the  range  Ms  =  1  to  10  and  =  ]<r*  to  1000  cm  Hg. 

The  data  in  the  present  report  cover  the  same  range  of  conditions 
in  imperfect  nitrogen.  The  data  of  Hilsenrath  and  Klein  (Ref.  4),  Little 


*The  following  terminology  will  be  used  in  this  report.  An  ideal 
gas  will  denote  a  gas  obeying  the  equations  p  -  pRT  and  h  =  CPT.  A 
perfect  gas  will  denote  one  obeying  p  =  ZpRT  which  includes  dissocia¬ 
tion  and  ionization  neglecting  intermolecular  effects.  An  imperfect  gas 
obeys  p  =  ZpRT  but  includes  intermolecular  forces.  The  term  real  gas 
is  reserved  for  a  gas  under  actual  (experimental)  conditions  and  is  not 
subject  to  the  theoretical  approximations. 
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and  Neel  (Ref.  21),  and  Grabau  and  Brahinsky  (Ref.  22)  were  interpo¬ 
lated  at  constant  pressure  and  entropy,  and  the  specific  heat  and  speed 
of  sound  data  were  computed  by  Neel  and  Brahinsky  for  the  single  set 
of  data  from  the  above  three  sources.  The  present  results  supplement 
the  earlier  calculations  of  Lewis  and  Burgess  in  the  range  Ms  =  1  to  6  . 
In  the  region  where  the  present  data  overlap  the  previously  published 
data  at  lower  pressures,  the  results  are  identical.  However,  the  data 
presented  herein  are  plotted  to  a  larger  scale  than  used  previously  to 
permit  increased  accuracy  in  reading  and  interpolating  the  graphical 
results. 

The  purpose  of  the  present  report  is  to  present  the  solutions  of  the 
normal  shock  wave  equations  where  the  effects  of  intermolecular  forces 
have  been  taken  into  account.  The  conditions  behind  the  incident  and 
reflected  shocks,  the  stagnation  conditions  upstream  and  downstream, 
and  conditions  immediately  downstream  of  a  standing  normal  shock 
wave  (see  Fig.  1)  are  presented  in  graphical  form.  Charts  for  pres¬ 
sure,  p  ,  temperature,  T  ,  velocity,  u  ,  enthalpy,  h  ,  compressibility 
factor,  Z,  entropy,  S,  and  speed  of  sound,  a,  in  the  above-mentioned 
regions  will  be  given  non-dimensionalized  by  appropriate  quantities  in 
the  quiescent  gas.  In  addition,  Ma  and  (q  \fTT ),0 '  >  are  also  given. 

SECTION  II 

NORMAL  SHOCK  EQUATIONS  AND  CALCULATION  PROCEDURE 

2.1  INCIDENT  SHOCK 

Consider  first  a  plane  shock  wave  moving  into  a  quiescent  gas  as 
shown  in  the  following  sketch: 


© 

»»  u2 

►  us 
© 

U1  =  0V— ' 

( 

a) 

The  shock  wave  has  a  velocity  Ug  with  respect  to  the  coordinate  system 
fixed  in  the  quiescent  gas.  By  superposing  a  velocity  (-Us)  on  the  co¬ 
ordinate  system,  the  shock  now  is  stationary  in  the  new  coordinate  sys- 

tern  as  shown  below: 

© 

O 

US  -  u2 

us 

(b) 


where  all  other  quantities  are  invariant  under  the  transformation. 
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The  continuity,  momentum,  and  energy  equations  for  the  system 
shown  in  sketch  (b)  are 


P2*(US  -  u2)  =  p*  Us 

(1) 

+  p*2  (us  -  02  =  P*  +  P*  usj 

(2) 

h2  l-  ‘4  (US  -  u,)J  =  h,  +  %  Us" 

(3) 

Eliminating  u2 

from  Eqs.  (1)  through  (3)  yields 

P3/p1  =  1  +  PiUs*  U  -  Pjjl/PjRT,, 

(4) 

h2/h,  =  1  +  USJ  (1  -  Pi2  )/2h, 

(5) 

Equations  (4), 

(5),  and  the  equation  of  state, 

P2  =  p(p» 

(6) 

form  the  system  of  equations  for  the  unknowns  p2 ,  p2 ,  and  h2  to  be  com¬ 
puted  with  the  given  shock  velocity.  Us,  and  the  conditions  in  the 
quiescent  gas  denoted  by  subscript  1.  Double  subscripts  denote  ratios, 

e-  &*  *  Pl2  55  Pi /P2 « 

The  thermodynamic  data  were  interpolated  at  constant  pressure  by 
Neel,  and  the  data  were  arranged  in  tables  of  the  form  <p  =  ( log  p ;  log  p , 

h/RT,  S/R ,  T,  Z)*  .  The  data  were  interpolated  in  the  range 
log  p  =  -7. 4(0. 1)4.8. 


The  calculating  procedure  for  the  solution  of  Eqs.  (4)  through  (6) 
was  as  follows:  With  the  given  conditions  and  an  initial  guess  for 

Pl}  (say  pffl,  Eqs.  (4)  and  (5)  were  solved  for  p2  and  h2 .  The  constant 
pressure  interpolated  thermodynamic  data  were  used  to  determine  p2 

by  Eq.  (6),  and  then  P<1  was  computed.  A  comparison  between  p  J1' 

12 

and  p  <0  was  made,  and  the  procedure  was  repeated  until 

“  f!  ^  I  <  0.00005  .  All  properties  except  the  speed  of  sound,  a, 

were  determined  from  a  table-look-up  procedure  in  the  constant  pres¬ 
sure  interpolated  data.  The  speed  of  sound  data  were  tabulated  in  the 
form  a  =  a  (S/R;  log  p,  log  p ,  h/RT ,  T ,  Z '  and  in  the  range  S/R  -  11.6(0.2)50. 


*This  notation  is  used  to  denote  log  p  and  any  one  of  the  remain¬ 
ing  quantities  as  independent  variables;  thus  any  quantity  shown  can  be 
determined  from  the  two  independent  variables. 
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2.2  STANDING  SHOCK 


A  standing  normal  shock,  such  as  the  bow  shock  in  front  of  a  body, 
is  shown  in  the  following  sketch: 


Uc 


© 


where  2s  denotes  conditions  immediately  downstream  of  the  standing 
shock. 


The  conservation  equations  are 


* 

*  ,  *  2 
P2  + 

h,  +  V2  UjJ 

and  eliminating  ujS  yields 


P  U 

r  2S  2S 

*  *  i 

P«  +  Pis  uis 

&S  +  % 


pJS/p2 "  1  +  C1  -  W<,«*)]/p.RT« 


his/lli  =  1  +  uaJ[l  -  (pt/p2S)]/ 2h} 

Equations  (10)  and  (11)  with  the  equation  of  state, 

P2S  =  p  (  p  ,  h  ) 


(7) 

(8) 

(9) 

(10) 

(ID 


(12) 


complete  the  system  of  equations  in  the  unknowns  p3s,  hjs  ,  and  pJS  •  The 
method  of  solution  was  identical  to  that  used  in  the  incident  shock  case. 


2.3  R  EFLECTED  SHOCK 

A  normal  shock  wave  reflected  from  the  closed  end  of  a  (one¬ 
dimensional)  tube  is  shown  in  the  following  sketch: 


which,  when  transformed  to  a  coordinate  system  fixed  with  respect  to 
the  shock  wave,  becomes 
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® 

UR  +  u2 


The  conservation  equations  in  the  transformed  coordinate  system  are 


P2  (Ur  u2)  =  p*Ur  ( 13) 

P*  +  P*  (Ur  +  u2)  =  p*  +  p*  Ur1  (14) 

h*  -i-  %  (Ur  +  =  h,  +  %  Ur2  (15) 

and  eliminating  u2  results  in 

Ps/Pa  -  1  +  P2 U r2  (1  -  p2S  )/plPJ  RTa  (16) 

h,/h2  =  1  +  PjiPjj  -  DU  +  p„  )  RTa/2p2h2  (17) 

with  the  equation  of  state  in  the  form 

ps  =  p(p , h )  ( 18) 

Using  Eqs.  (1)  and  (13)  to  eliminate  u2  gives  an  expression  for  UR  in 
terms  of  Ug  ,  pl2  ,  and  p25  as  follow's : 

Ur  -  US(1  -  pu)p„/i\  ~  P„)  (19) 


Equations  (16)  through  (19)  were  solved  in  a  similar  manner  to  the 
incident  shock  case. 


2.4  STAGNATION  CONDITIONS 

The  isentropic  stagnation  conditions  upstream  and  downstream  of  a 
standing  shock  wave  were  determined,  respectively,  from  table-look-up 
procedures  of  the  form  4>  =  <4(Ss,  h0)  and  6  =  (S2S,  h0)  where  <b  and  ^ 

represent  ail  thermodynamic  quantities  and  hD  =  h2  +  |/j  u/  . 

SECTION  III 

RESULTS  AND  DISCUSSION 

The  results  of  the  calculations  are  presented  in  Figs.  2  through  6.  * 
The  ideal  gas  (y  =  1.4)  results  are  also  shown  for  comparison.  Figure  2 


*The  tabulated  numerical  results  can  be  obtained  from  the  first 
author  upon  written  request. 
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shows  the  pressure,  density,  temperature,  enthalpy,  velocity,  speed 
of  sound,  Mach  number,  compressibility  factor,  and  entropy  behind  the 
incident  shock  wave.  Similar  data  for  the  conditions  downstream  of  the 
standing  shock,  stagnation  conditions  upstream  of  the  shock,  and  con¬ 
ditions  behind  the  reflected  shock  are  given  in  Figs.  3  through  6. 


The  sphere  stagnation  heat-transfer  rate  (q  v' rn  )20,  was  computed  by 
the  Fay  and  Riddell  (Ref.  24)  formula  similar  to  the  form  used  previ¬ 
ously  by  Lewis  and  Burgess  (Refs.  25  and  26) 


qyrn  =  7.55274  x  10“  *  [pw  pw  )°*‘  (pjg  'P2B  -)*'*  (hc  -  Kw  )  [(  pj0  , /pM  ,)  (1  -  P,/pJcj] 


with  dimensions  Bui  V  in. /ft1 -sec  .  The  viscosity  data  were  those  of  Ahtye 
and  Peng  (Ref.  27). 


Where  the  thermodynamic  properties  could  be  interpolated  from 
the  data  of  Hilsenrath  and  Klein  (Ref.  4),  and  the  data  of  Grabau  and 
Brahinsky  (Ref.  22),  conditions  were  computed  for  incident  shocks  in 
the  range  Ms  =  1  (0.5)  JO  and  quiescent  gas  pressures  in  the  range 
P)  =  1.32  x  10-6  to  13.16  atm  . 

The  quantities  most  affected  by  the  intermolecular  forces  are 
Z ,  p ,  and  a  .  The  effects  on  Z  above  one  atmosphere  in  region  2  can 
be  seen,  for  example,  from  Fig.  2h.  The  effect  on  p  was  demonstrated 
earlier  by  Lewis  and  Burgess  (Ref.  15)  and  Lewis  (Ref.  17)  by  compari¬ 
son  of  their  results  with  the  perfect  gas  results  of  Bernstein  (Ref.  19) 
and  Ahtye  and  Peng  (Ref.  27).  The  effect  on  a  was  also  demonstrated  by 
Lewis  and  Neel  (Ref.  13)  by  similar  comparisons  with  the  perfect  gas 
speed  of  sound  data  of  Ahtye  and  Peng  (Ref.  27). 

Comparisons  were  previously  made  by  Lewis  and  Burgess  (Ref.  15) 
between  their  results  and  those  of  Bernstein  (Ref.  19)  at  Fi  up  to 
2fif)  cm  Hg  .  Below  one  atmosphere,  the  differences  were  well  within 
plotting  accuracy,  i.  e. ,  less  than  one  percent.  At  2.  6  atmospheres  in 
the  quiescent  gas,  the  differences  in  Z5  ,  ps  ,  and  Ur  were  as  much  as 
10  percent  at  Y1S  =  12  .  The  influence  of  the  intermolecular  forces  was, 
of  course,  largest  at  the  highest  densities.  Lack  of  perfect  gas  results 
at  higher  pressures  prevented  further  comparisons,  and  similarly, 
lack  of  perfect  gas  speed  of  sound  data  and  results  in  the  20  region  pre¬ 
vented  further  comparisons. 


The  results  presented  herein  at  Pi  =  1000  cm  Hg  were  almost  entirely 
based  on  the  data  of  Grabau  and  Brahinsky  (Ref.  22),  which  were  ob¬ 
tained  by  extrapolating  published  thermodynamic  data  to  1000  amagats 
and  by  demanding  internal  consistency  and  satisfaction  of  the  laws  of 


6 


AF  DC-T  R-66-1 16 


thermodynamics.  However,  the  data  are  subject  to  approximate 
numerical  treatment,  and  thus  the  results  contained  herein  reflect  the 
inaccuracies  caused  by  numerically  fitting,  extrapolating,  interpolating, 
and  differentiating  the  basic  thermodynamic  data. 

The  effects  of  joining  the  extrapolated  data  of  Grabau  and  Brahinsky 
with  the  previous  data  of  Hilsenrath  and  Klein  (Ref.  4}  can  be  seen  in 
Fig.  5f.  The  curves  at  200  and  loo  cm  Hg  below  Ms  =  7.5  were  calculated 
based  on  the  extrapolated  data  of  Grabau  and  Brahinsky,  whereas  above 
Me;  =  7.5  the  data  are  those  of  Hilsenrath  and  Klein.  The  differences 
where  the  data  sources  join  are  small,  and  the  trends  of  the  results  are 
as  expected.  Since  Hilsenrath  and  Klein  included  only  two  virial  coeffi¬ 
cients,  whereas  the  extrapolation  of  Grabau  and  Brahinsky  was  not 
based  on  a  virial  equation  of  state,  the  extrapolated  Z  must  be  equal  to 
or  greater  than  the  one  including  only  the  second  virial  correction. 

Small  deviations  are  thus  to  be  expected  in  the  present  results,  but  it 
is  hoped  that  the  results  are  little  affected  by  the  numerical  procedures. 

Finally,  comparison  of  the  results  in  this  report  with  those  in  air 
(Ref.  23)  for  the  same  initial  conditions  would  indicate  the  effects  of 
oxygen  vibration  and  dissociation  on  the  gasdynamic  properties  behind 
incident  and  reflected  shock  waves.  The  effects  are  larger  at  the  lower 
P,  where  the  effects  of  oxygen  dissociation  can  easily  be  seen  from 
Ref.  23  for  quantities  such  as  p,  T,  and  a  in  the  stagnation  and  reflected 
regions.  It  is  thus  hoped  that  the  present  charts  will  be  useful  to  those 
interested  in  low  as  well  as  high  initial  pressure  experimental  and  theo¬ 
retical  studies  in  nitrogen. 


SECTION  IV 
CONCLUDING  REMARKS 


Normal  shock  wave  parameters  have  been  calculated  using  im¬ 
perfect  gas  thermodynamic  properties.  The  results  well  below  one 
atmosphere  pressure  in  the  quiescent  gas  are  in  agreement  with  the 
previous  calculations  using  pei'fect  gas  properties.  At  one  atmosphere 
and  above,  the  effects  of  the  intermolecular  forces  on  the  imperfect  gas 
properties  and  hence  on  the  normal  shock  parameters  were  previously 
demonstrated  (see  Ref.  15).  Moreover,  as  increases  imperfect  gas 
effects  increase. 
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